home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / spptrs.z / spptrs
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          SSSSPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPPTRS - solve a system of linear equations A*X = B with a symmetric
  10.      positive definite matrix A in packed storage using the Cholesky
  11.      factorization A = U**T*U or A = L*L**T computed by SPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SPPTRS( UPLO, N, NRHS, AP, B, LDB, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDB, N, NRHS
  19.  
  20.          REAL           AP( * ), B( LDB, * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      SPPTRS solves a system of linear equations A*X = B with a symmetric
  24.      positive definite matrix A in packed storage using the Cholesky
  25.      factorization A = U**T*U or A = L*L**T computed by SPPTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangle of A is stored;
  31.              = 'L':  Lower triangle of A is stored.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      NRHS    (input) INTEGER
  37.              The number of right hand sides, i.e., the number of columns of
  38.              the matrix B.  NRHS >= 0.
  39.  
  40.      AP      (input) REAL array, dimension (N*(N+1)/2)
  41.              The triangular factor U or L from the Cholesky factorization A =
  42.              U**T*U or A = L*L**T, packed columnwise in a linear array.  The
  43.              j-th column of U or L is stored in the array AP as follows:  if
  44.              UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j; if UPLO =
  45.              'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
  46.  
  47.      B       (input/output) REAL array, dimension (LDB,NRHS)
  48.              On entry, the right hand side matrix B.  On exit, the solution
  49.              matrix X.
  50.  
  51.      LDB     (input) INTEGER
  52.              The leading dimension of the array B.  LDB >= max(1,N).
  53.  
  54.      INFO    (output) INTEGER
  55.              = 0:  successful exit
  56.              < 0:  if INFO = -i, the i-th argument had an illegal value
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.